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The uncertainty relation for the photons in three dimensions that overcomes the difficulties caused 
by the nonexistence of the photon position operator is derived in quantum electrodynamics. The 
photon energy density plays the role of the probability density in configuration space. It is shown 
that the measure of the spatial extension based on the energy distribution in space leads to an 
inequality that is a natural counterpart of the standard Heisenberg relation. The equation satisfied 
by the photon wave function in momentum space which saturates the uncertainty relations has 
the form of the Schrodinger equation in coordinate space in the presence of electric and magnetic 
charges. 



The uncertainty relations that first appeared in the 
fundamental Heisenberg paper on the conceptual con- 
tent of quantum mechanics [l| play an important role 
in elucidating specific properties of the quantum world. 
The aim of this Letter is to derive in the framework of 
quantum electrodynamics a sharp uncertainty relation 
for photons treated as three-dimensional objects. The 
photon wave function that saturates our inequality has 
already appeared in different contexts. 

An uncertainty relation for the momentum and the 
center of mass coordinate along a given direction can 
easily be obtained from the commutators between the 
generators of the Poincare group. It follows from general 
principles Q that the commutators between the compo- 
nents of the momentum and the Lorentz boost are: 



ihSijH. 



(1) 



Therefore, each component of the center of mass coordi- 
nate (the boost operator divided by the energy operator 
with proper symmetrization to secure Hermiticity), 



Ri 



(2) 



is a canonically conjugate variable to the corresponding 



. Hence, the 
/A follows 



component of momentum, [Ri,Pj\ — ih5, 
standard Heisenberg relation AxjApf > 
separately for every direction, with Gaussian functions 
saturating each inequality. These simple one-dimensional 
uncertainty relations hold for every relativistic system. In 
the special case of photons they were obtained by Holevo 
in the framework of estimation theory. Holevo also 
proved that "the transversality condition precludes the 
three inequalities from becoming equalities simultane- 
ously." This was to be expected because the existence 
of one wave function that saturates all three standard 
Heisenberg uncertainty relations would mean that pho- 
tons are not any different from nonrelativistic particles. 
In a general approach based on the properties of the 



Poincare group, the impossibility to saturate simulta- 
neously the three uncertainty relations follows from the 
noncommutativity of the components of R, 



AUi - R,Pi + RiP, H 



(3) 



The expression in parentheses on the right side represents 
the spin: the total angular momentum minus the orbital 
part. Schwinger [43 has shown that for massless particles 
with helicity A the right-hand side in (jS]) becomes: 



Rx R 



-iXk/k^ 



(4) 



He also found a very rough estimate of the lower bound 
in the three-dimensional uncertainty relation. 

In what follows we fully realize the Schwinger program 
by deriving a precise form of the uncertainty relation for 
photons in three dimensions. We also find the photon 
wave function in momentum space that saturates the in- 
equality and we confirm the mathematical analogy with 
magnetic monopoles discovered by Schwinger. Our un- 
certainty relation for photons has the form: 



ArAp > 4h, 



(5) 



and is saturated (i.e. becomes an equality) by — what 
might be called — the fundamental solution of Maxwell 
equations. The lower bound in this relation is 4 times 
larger than the Schwinger estimate. It is also larger than 
3/2h that would be obtained by blindly extending the 
one-dimensional relation for photons to three dimensions. 
The larger value is due, of course, to the impossibility to 
saturate simultaneously the one-dimensional uncertainty 
relations for photons in all three directions. 

In order to overcome the problem of the nonexistence 
of a bona fide photon position operator and the associ- 
ated probability density in coordinate space we rely on 
the photon energy density. To measure the photon spread 
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in coordinate space we use the second moment of the pho- 
ton energy distribution — the generator Kq of a conformal 
transformation [Hj. 

There is no problem with the momentum distribu- 
tion of photons since there exists a well defined photon 
wave function in momentum space. This function ap- 
pears in the representation theory of the Poincare group 
expounded by Wigner @. In order to describe com- 
pletely the state of a photon we need two such functions, 
say /±(fc), that refer to two helicities (two circular po- 
larizations). Under all proper Poincare transformations 
these two functions transform separately but they are in- 
terchanged under reflections. In quantum theory these 
functions (normalized to 1) become the probability am- 
plitudes for finding the photon with momentum k and 
the helicity A = +lorA = — 1. Following Wigner, we 
choose the relativistic form of the norm, 



(6) 



where we allow for both polarizations to be present. 

The straightforward Fourier transformation of the pho- 
ton wave function in momentum space does not produce 
a local function in coordinate space 0, S] ■ To describe 
the localization properties in coordinate space we shall 
rely solely on the electromagnetic fields. 

It is convenient to express the energy density in terms 
of a complex combination of the electric and magnetic 
field operators, named in Q the Riemann-Silberstein vec- 
tor. 



F{r,t) 



D{r,t) ^■Bir,t) 



2fi 



(7) 



The energy density operator is equal to F^{r,t)-F{r,t). 
We shall take the following expectation value: 



he 



(frr^ : F\r,t)- F{r,t) :) 



(8) 



as a measure of the spread in the coordinate space. In 
this way we follow closely the concept of photon local- 
ization underlying the theory of photodetection lol 
According to this theory, the localization of photons is 
determined by the correlation functions of the field oper- 
ators. Normal ordering in (|8]) eliminates the contribution 
from the vacuum state. We have chosen the origin of the 
coordinate system at the center of mass ((r) = 0). De- 
spite the presence of , Ar has the dimension of length 
like its counterpart in the Heisenberg uncertainty rela- 
tion. 

In quantum electrodynamics F{r,t) has the following 
representation in terms of the annihilation and creation 
operators for photons with positive and negative helici- 



ties a+(fe) and al(fc) (cf. for example, 0, Q), 
F{r,t) 



(27r)3/2 

X e(fe) a+{k)e'''-'-"^* + al{k)e~'''-''+"^* 



(9) 



The normalized polarization vector e{k) is taken as: 

e{k) = ^ I -kyk, - ikK I . (10) 



kl 



The polarization vector e{k) always distinguishes a cer- 
tain direction. We have chosen here this direction along 
the z axis. The commutation relations between the cre- 
ation and annihilation operators. 



ax{k),a{,{k')\ = 5xx'k6^^\k - k'), (11) 



follow from the commutation relations for field operators. 

We shall evaluate Ar in the one-photon state that is 
created from the vacuum state by the following combi- 
nation of creation operators: 



!./> = 



(Pk 



U{k)al{k) + f-{k)al{k)\ |0). (12) 



The norm of this state vector is given by Eq. ([6]) . We will 
not assume now that the wave functions are normalized 
to 1 since they will serve later as variational variables. 
Our measure of the uncertainty in position for this state 
reduces to the second moment of the classical energy den- 
sity (divided by he), 



^""/MW jd'rr'F*[r,t)-F{r,tl 



(13) 



where the c-number Riemann-Silberstein vector F{r,t) 
is built from the photon wave functions /±(fe), 



F{r,t) 



dPk 



(27r)3/2 

X e(fe) [/+(fe)e''^"-'"* + r_{k)e-'''-''+"^'] . (14) 

We may choose i = 0, since there is no preferred origin 
of time. Then, Ar can be converted into an integral in 
momentum space. 



-|/a 



(15) 



where D\ is the covariant derivative (in momentum 
space) on the light cone [1, [H, [l3| . The helicity A plays 
here the role of the charge, 

Dx = V -i\a.{k), (16) 
^(fe)- fc(fc2^;fc2) {-kyM). (17) 
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In the derivation of (|15l) we used the formula: 
, , n- f d^k 

X e^'=-a. (e,(fc)/+(fc)) - e-''="-9. (e,(fc)r(fc)) ] , (18) 
and the following properties of the polarization vector: 



ie* {k) ■ die{k) = ai{k) 
d,e*{k)-d,e{k)= ^ 



jU2 I jU2 " 



e{k)-d.,e{k) ^{), (19) 
a,e(fe)-a,e(fe) = 0. (20) 



The vector a{k) plays the role of the connection on 
the light cone that defines the parallel transport for the 
photon wave function in momentum space. When the 
phase of f\{k) is changed by exp(iA(/ii(fc)), the phase 
of e{k) must undergo a compensating change, e(fc) — >■ 
exp[— i(/!)(fc)]e(fe), to keep the field F{r,t) intact. This 
results in a gauge transformation (y.{k) — > Q;(fe)+V0(fc). 
Note, that the covariant derivative vector iD\ obeys the 
same commutation relations as the operator R, 



The existence of a finite lowest value of the product 
ArAp is an expression of the uncertainty relation for the 
photon. This value will be calculated by a variational 
method. 

We test this method first for the standard Heisenberg 
uncertainty relation. Again, we choose the center of mass 
and the center of momentum frame, so that (r) = and 
(p) = 0- The expression to be minimized is: 



9 9 
2 '^r<^p 

W = 



h? Jd-^r ip*{r)ip{r) Jd-^r tlj*{x)ip{x) 

(24) 

Variation with respect to ip*{r) leads to the Schrodinger 
equation for the three-dimensional harmonic oscillator, 



i [-crlA + r^al] V'(r) = ^|J{r) 



(25) 



where we assumed at the end that the wave function is 
normalized to 1. The operator appearing on the left-hand 
side has the eigenvalues zu{n+5/2). Comparing this with 
the right-hand side, we obtain zu = (n + 3/2). The lowest 
value of t37 is 3 /2 and it leads to the standard Heisenberg 
uncertainty relation in three dimensions. 



iDx X iDx = iW x a{k) = -iXk/k^ 



(21) 



Therefore it can be viewed as a representation of R in 
momentum space. 

The last equation means that a(fc) can be treated as a 
vector potential of a magnetic monopole. Such a vector 
is always singular on the Dirac string [31 • Our choice of 
the polarization vector e(fc) leads to the string along the 
whole z axis {cf. for example, Eq. (16) of Ref. ^ISi]}. 

Our uncertainty relation for photons involves the prod- 
uct ArAp like the standard Heisenberg relation. The 
natural choice for Ap is: 



Ap = 



■E 



h\k-{k)\ \h{k)\' 



(22) 



To simplify the calculations we shall choose the coordi- 
nate system comoving with the wave packet in such a 
way that h{k) = (center of momentum frame). Then, 



Ap : 



■E 



d^k 



k\h{k)\^ ^h{k). (23) 



(26) 



The wave function of the ground state — the Gaussian — 
saturates this inequality. Owing to the linearity of the 
variational equation (|25p we did not have to use a La- 
grange multiplier. Note that for the wave functions that 
saturate the inequality the uncertainties in position and 
momentum are evenly distributed, when measured in 
conjugate units, namely cr^ = \/3/2 a, ap = ■\/3/2 h/a. 

Now, we apply an analogous variational procedure to 
the dimensionless expression 7 = ArAp/h. With the use 
of ([T6| and (fT7|) we rewrite the formulas (fT5|) and (f23|) in 
spherical coordinates. 



Ar : 



1 



J2 dkk^ dOsine dip \dkf. 



fc2 sin^ 



^P=TrTTHy2 dkk^ d9 sine d^\f, 
11/11 , Jo Jq Jq 



\deh\^ 
k^ 

(27) 
(28) 



Variation of 7 with respect to fx{k,9,(p) produces the 
following equation: 



1 



-dgsin9dg 



1 + 2iX cos 



27 



+ 7 



(29) 



where we assumed (after performing the variation) that the norm of / is equal to 1 and we introduced a dimen- 
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sionless variable k = hk/ /^p. After the separation of 
variables, f\{K.,9,Lp) — IC{k)Q\{9) exp^irrup)^ we obtain: 



1 



27 

K 



IC(k) = -7/C(k), (30) 



sm ( 



■ de sin 9 da 



vn? + }? ~ 2\m cos 9 



sm 



em 



^j{3 + l)Qx{e). (31) 

The equation for Q\{9) is exactly the same as in the 
theory of magnetic monopoles (cf. for example, [islflGj). 
We wrote instead of 1 to have the form of Eq. (PT|) 
identical with the monopole case. In this equation A 
plays the role of the product of the electric and magnetic 
charges divided by he. The presence of the monopole 
string along the z axis excludes the value j = Q. The 
contributions to Ar and to Ap from the wave functions 
with A = ±1 are equal and independent. For definiteness, 
we choose A = 1. 

The equation for the radial part has the same form 
as the Schrodinger equation for the Coulomb problem. 
Hence, the spectrum of the operator appearing on the 
left-hand side of Eq. ^ is -7^/(nr + j + 1)^. Since this 
must be equal to —7, we obtain 7 = (n^ +j + 1)^. We are 
interested in the smallest value of 7 which corresponds 
to rir = 0, j = 1. This gives 7 = 4 as in Eq. ([5]). 

There are three independent solutions of Eq. (|3T1) for 
J = 1 that saturate the uncertainty relation. In the angu- 
lar momentum basis they correspond to m = 0, ±1 but 
we found the Cartesian basis more convenient. Normal- 
ized solutions of the variational equation that saturate 
the uncertainty relation are: 



f^^'\k,9,ip)^^ksm9e-'"', (32a) 
f^''\k,9,ip) = -^k(-ism(p-cos9cos(p)e~''"', (32b) 
f^y\k,9,ip) = — =k{i cos LP -cos 9 sin Lp)e-'"', (32c) 



where a is an arbitrary scaling parameter. These three 
solutions are not essentially different. They correspond 
to three choices of the polarization vector e{k) obtained 
from (ITUl) by cyclic permutations of the coordinate axes: 
kz ^ kx ^ ky V kz • 

The uncertainty relation ([5]) does not depend on a. 
Moreover, the uncertainties in position and momentum 
for the wave functions that saturate the inequality are 
again, as is the case for the Heisenberg relation, equally 
distributed: A,. = 2a, Ap = 2h/a. The larger value of 
the numerical coefficient in the uncertainty relation for 
photons than in the Heisenberg relation ((26|) is due to 
the magnetic monopole term in Eq. ([3T|). Since this term 
excludes s waves, the maxima of the functions p2p are 
shifted away from zero. 



The solution of Maxwell equations that saturates the 
uncertainty relation deserves to be named the fundamen- 
tal solution. It can be built from the second derivatives 
(Whittaker construction [l3l ) of the complexified funda- 
mental solution $ of the d'Alembert equation. 



<i>(r, r) = [r^ - (t - iaf] 



(33) 



where t = ct. For our choice of e(fc), the wave function 
/(^) leads to: 



F(^)(r,T) 
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47r2 



dxdz ^ 
dydz - 



idydr 



Uy 



(34) 



The superscript (z) indicates that the z axis was sin- 
gled out in the definition of the polarization vector. The 
remaining two solutions, built from Z^^-* and /^^\ are ob- 
tained from p4|) by cyclic permutations of partial deriva- 
tives. We may also apply dual rotations to further en- 
large the set of solutions of Maxwell equations that sat- 
urate the uncertainty relation. This whole set may be 
viewed as just one fundamental solution acted upon with 
various symmetry transformations. At the moment of 
maximal compression {t — 0) they all have the same 
spherically symmetric energy density £: 



he 



7r2 (a2+r2)4 



(35) 



The simplest solution, given by the expression (|34|) . re- 
duces at i = to a pure electric field: 



4a2 {2zx — 2ay, 2yz + 2ax, 



(36) 



whose all field lines are ellipses. 

The fundamental solution of Maxwell equations ap- 
peared in various investigations. Its exceptional prop- 
erties were recognized by Synge who called it "an 
electromagnetic model of a material particle." Later, this 
solution appeared as the main ingredient in the construc- 
tion of electromagnetic knots by Hopf fibrations 



ia-2 



and in the study of vortex lines embedded in the solutions 
of Maxwell equations [i^. It was even used to explain 
the ball lightning [24 1 . 
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